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The stochastic series expansion quantum Monte Carlo method is used to study thin ferromag- 
netic films, described by a Heisenberg model including local anisotropies. The magnetization curve 
is calculated, and the results compared to Schwinger boson and many-body Green's function cal- 
culations. A transverse field is introduced in order to study the reorientation effect, in which the 
magnetization changes from out-of-plane to in-plane. Since the approximate theoretical approaches 
above differ significantly from each other, and the Monte Carlo method is free of systematic errors, 
the calculation provides an unbiased check of the approximate treatments. By studying quantum 
spin models with local anisotropies, varying spin size, and a transverse field, we also demonstrate 
the general applicability of the recent cluster-loop formulation of the stochastic series expansion 
quantum Monte Carlo method. 

PACS numbers: 75.10. Jm,75.40.Mg,75. 70.Ak,75. 30. Gw 
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I. INTRODUCTION 

The driving force behind much of the research on thin 
magnetic films is their application in data storage de- 
vices. Magnetic thin films also display many remarkable 
physical phenomena, such as the reorientation effect ,EJ in 
which the axis of magnetization changes as a function of 
film thickness, temperature and applied fields. Many the- 
oretical methods, such as ab-initio calculations^ mea 
field theories,Ba classical Mflftf^r^Pl^'-' simulations! 
Green's-, function methodsQEHall3llJ't3 and Schwinger 
bosongll3 have been applied to thin- film systems. Much of 
this work has focused on using the Heisenberg model to 
study the reorientation effect. The ground state and the 
lowest (one-magnon) excitations are known for the two- 
dimensional ferromagnetic Heisenberg model, but there 
is no closed form analytic solution at finite temperatures. 
Since results obtained using different approximate meth- 
ods differ a great deal from each other, there is a need 
for an unbiased check of the various methods used; this is 
achieved by quantum Monte Carlo (QMC) calculations. 
In Ref. ^, for exarpfjle. it has been shown by comparing 
with QMC resultfJijIla that the TyablicovB decoupling 
(random phase approximation: RPA) is a very good ap- 
proximation for the magnetization of a spin 5 = 1/2 
monolayer in an external magnetic field (perpendicular to 
the film plane) . The main purpose of the present paper 
is to show the feasibility of large-scale quantum Monte 
Carlo (QMC) calculations, free of systematic errors, in 
the study of thin magnetic films. In order to achieve this 
goal we have included higher spins, local anisotropies, 
and a transverse magnetic field in the operator-loop for- 
mulation-of the stochastic series expansion (SSE) QMC 
method.EII 

Two often-cited trends in quantum Monte Carlo de- 
velopment are the emergence of methods free of system- 



atic errors 



and thi 



cient loop-cluster algorithms 



of highly efR- 
In this work 



we particularly want to emphasize the general applica- 
bility of the SSE operator-loop method, which makes it 
possible to use the same algorithm to study a wide variety 
of Hamiltonians. Whereas previously it was necessary to 
rewrite large sections of the computer code when chang- 
ing the model, one can now use the same code (compiled 
only once) to simulate a wide range of different systems. 
The user of the program no longer necessarily needs de- 
tailed knowledge of the algorithm and code to be able to 
conduct a thorough study of many quantum spin models 
in any dimension with, for example, non-zero magnetic 
fields, anisotropies and varying spin size. 

A brief description of the QMC method is given in 
Sec. II. Thereafter we discuss the applicability of the 
method and the introduction of general spin size and a 
transverse field in the SSE operator-loop algorithm. In 
Sec. [H we compare some examples of the QMC simula- 
tions to approximate theoretical approaches. Finally, we 
comment on possible future applications of QMC in the 
context of thin magnetic films. 



II. SSE CLUSTER-LOOP ALGORITHM 

There ace— jexcellent descriptions of the SSE loop 
algorithm,EZrE3 so we only give a brief summary here in 
order to introduce the general framework of the method. 
We will, however, try to describe the main features of the 
method pictorially. The focus is on the new aspects that 
arise when introducing arbitrary spin size and a trans- 
verse field. 

We consider a lattice spin model described by a Hamil- 
tonian H. The SSE method relies on a Taylor expansion 
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of the partition function Z: 
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where |a) are basis states in which the matrix element 
above can be evaluated, and /? is the inverse temperature. 

To describe the updating procedure, we write the 
Hamiltonian as a sum over all M bonds representing in- 
teracting spins in the system 



M 
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The bond operator Hi, can be decomposed into its diag- 
onal and ofF-diagonal parts. 
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where subscript D denotes a diagonal operator and O an 
off-diagonal operator. For a ferromagnetic Heisenberg 
model these two operators are of the form 



and 



HD,b — S'j^(;,)5'J(f,) 



Ho.b - St{b)^j{b) + ^t{b)^t{b)^ 
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where i{b) and j{b) denote the two spins connected by 
bond b. If we introduce a cutoff order L in the Taylor 
expansion (which, when done properly, does not cause 
any systematic errorsEll) , and include additional unit op- 
erators /, the expansion can be rewritten in the form 



H = 



EE 



\SL\a), 
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where Sl is an operator string 
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with Hp e {H]:) ii,Ho,b,I}- The Monte Carlo procedure 
must sample the space of all states \a) and all operator 
sequences Sl- 

We next consider the SSE space in more detail. De- 
noting a propagated state by 
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the matrix element in Eq. (^) can be written as a prod- 
uct of elements of the form {a{p)\Hi,\a{p — 1)), where the 
bond-operator Hh only acts on two spins. For the spin- 
1/2 Heisenberg model the only elements that can appear 
are (T,T \Ht\ T,T), (T, i \Hb\ ti), (T,i \Hb\ i,T), and 
the spin-reversed versions of the same set. From now on 
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FIG. 1: An SSE configuration for a four-site spin-1/2 model. 
A filled (empty) circle denotes a state with spin up (down). 
A horizontal bar indicates a vertex, corresponding to a bond 
operator Hb, labeled on the right. The propagated states 
\a{p)) are labeled on the left side. 
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FIG. 2: A possible outcome of performing a diagonal update 
on the configuration in Fig. |l| 



we refer to these matrix elements, consisting of four spin 
states and a bond operator, as "vertices". The matrix 
element in Eq. (^) can thus be viewed as a list of ver- 
tices of the above kind. In Fig. |l] we depict such a list 
graphically. 

In the operator-loop algorithm, two basic updates en- 
sure that the complete SSE space is sampled. The di- 
agonal update attempts to exchange diagonal operators 
HD,b and unit operators /. The probability for inserting 
a diagonal operator (exchanging it for a unit operator) 
at position p in the operator sequence is 



AIl3{a{p)\HD^b\a{p)) 
L — n 



(9) 



while the probability for removing a diagonal operator is 



L-n + 1 



M(3{a{p)\HD,bHp)) 



(10) 



In a diagonal update, one exchange attempt is made for 
each diagonal and unit operator. A typical outcome of a 
diagonal update is shown in Fig. ^ 

The second type of update is a global operator-loop 
update, which leaves unit operators unaffected. The idea 
of the loop move is to form and flip a closed loop of spins 
in the vertex list. In the process both the affected ver- 
tices and states are changed. The operator-loop update 
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(a) (b) 




(c) (d) 

FIG. 3: A loop update. Unit operators have been removed 
from the configuration in Fig. bl Vertical lines (links) show 
how the spins at difi^erent vertices are connected. Short hori- 
zontal bars denote link-discontinuities. 



together with the above diagonal update ensure that the 
complete SSE configuration space is sampled. 

The loop move is easy to formulate. In Fig. js] we illus- 
trate a loop move graphically. We have removed all the 
unit operators / from the operator sequence in Fig. 
Furthermore we only show the spins that are members of 
a vertex, and we turn the configuration into a linked list 
by connecting the same spin at different vertices by ver- 
tical lines, which we refer to as "links". Next we describe 
the loop move in more detail. 

First a random spiii, belonging to a vertex, is selected 
and flipped, see Fig. ||(a). We refer to this spin as the 
"entrance" spin to the first vertex. For the loop to pro- 
ceed, we select an "exit" spin, also belonging to the same 
vertex, see Fig. ^(b). The probability of choosing a given 
exit spin is proportional to the vertex matrix element, 
that results from flipping the entrance and exitJegs. This 
choice of exit spin ensures detailed balanceoThe exit 
spin is linked to another spin in the linked vertex list, 
and this spin is chosen as the entrance spin to the next 
vertex, see Fig. ^(c). In this manner a loop of spins on the 
space-time lattice is formed and flipped. Note that when 
the first spin is flipped a discontinuity is introduced into 
the vertex list. A link-discontinuity appears when two 
spins linked together are not in the same state. When, 
for example, the initial entrance and exit spins (assum- 
ing they are different) are flipped, two discontinuities are 
introduced in the linked list, see Fig. ^(b), where a hnk- 
discontinuity is marked by a short horizontal bar. One of 
the discontinuities gets propagated by the loop until the 
loop passes through the initial spin a second time, when 
the discontinuities "annihilate" each other and the loop 



closes, see Fig. ^(d). A practical and sufficient criterion 
for closing the loop is that the exit spin (after having 
been flipped) is in the same state as the next entrance 
spin (before flipping it). When this is the case, all discon- 
tinuities have been removed and the loop can be closed. 



A. Applicability of the algorithm 

The SSE loop algorithm is extremely general and can 
be efficiently applied to any model with a positive definite 
partition function (if this is not the case, one epcounters 
the usual difficulties due to the sign problemE3). This 
means that the algorithm can be applied to ferromagnetic 
and antiferromagnetic (on bipartite lattices) Heisenberg 
models of any quantum spin size in any dimension, in- 
cluding magnetic fields and local as well as exchange 
anisotropies. The interactions need not be short-ranged, 
as long as they do not cause frustration (which leads to 
the sign problem). 

From a programming point of view, one need only cal- 
culate the energies of all possible vertices (which typically 
is a simple task) and generate a table of exit probabilities 
as a function of four variables: the initial vertex type, the 
entrance spin, the entrance spin state and the exit spin. 

As an illustration we have written a code that works 
for a general d-dimensional Heisenberg model of the form 

i 

-J2[B^S! + B-Sf], (11) 

i 

where the exchange interaction is given by Jy (negative 
for antiferromagnets and positive for ferromagnets). S 
denotes the quantum spin and K2 and K4 are the sec- 
ond and fourth-order local uniaxial anisotropies. B is an 
external magnetic field. In this study we limit ourselves 
to a nearest neighbour interaction, but in principle one 
can include long-range interactions. In the following we 
describe the modifications to the loop algorithm that al- 
low the introduction of general spin size and a transverse 
field. 



B. Inclusion of arbitrary spin 

For the sake of simplicity, this-SSE loop-operator algo- 
rithm was originally describedlJ for a spin-1/2 system, 
but the generalization to a spin-S' system is straight- 
forward. For an arbitrary spin system there are more 
allowed vertices, but the probability for inserting and ex- 
tracting diagonal operators is still given by Eq. (|^) and 
Eq. (Iln) respectively. The only necessary change in the 
loop-update is the choice of the initial spin state. For 
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the spin-1/2 model the initial entrance spin was sim- 
ply flipped, but for the spin-S" model one can randomly 
choose among the 25* spin states that differ from the ini- 
tial spin state. All other aspects of the algorithm remain 
the same. 



C. Inclusion of a transverse field 

Previously a transverse field has-been included in the 
SSE algorithm using local updatesJlj In the present work 
we include the transverse field directly in the global 
loop update. There are very likely many different ways 
to formulate a loop-update which includes a transverse 
magnetic field. Here we will not make an exhaustive 
study, but rather present one possible algorithm, which 
we found particularly simple to implement. It is proba- 
bly not the most efficient update, and we are currently 
making a more detailed study of more efficient algorithms 
for including the transverse field. 

Including a transverse magnetic field is more involved 
because the total z component of the magnetization is 
no longer a conserved quantity, \H, ^ ■ Sl\ ^ 0, due 
to the presence of single lowering and raising operators 
(25'^ = + S~) in the Hamiltonian. As a consequence 
additional "flip vertices" , with different magnetization in 
the initial and final states (such as (t,T l-f^fcl Tii))? are 
also allowed in the vertex- list. The presence of flip ver- 
tices implies that, in general, one can exit the same exit 
spin with several different spin states. The increase in 
the number of exit states means that the exit probabil- 
ity is a function of five variables instead of four: original 
vertex, entrance spin, entrance state, exit spin and exit 
state. For a model that conserves the total z component 
of the magnetization the exit state is a function of the 
other four variables. 

To include a magnetic field, we also have to recon- 
sider how to start, and end, the operator-loop. With no 
magnetic field we introduce an initial link-discontinuity 
by flipping the first entrance spin. A second link- 
discontinuity appears when the first exit spin is flipped 
(unless the entrance and exit spins coincide, in which case 
the loop closes immediately). The second discontinuity 
is propagated as the loop progresses until the two discon- 
tinuities annihilate each other and the loop closes. With 
flip vertices allowed, the number of link-discontinuities 
can change by the introduction or removal of single flip 
vertices. If one therefore chooses to start, and end, the 
loop in the same way as previously described, the loop 
can close without having corrected for the initial discon- 
tinuity. In such a case it would be necessary to return 
to the starting point of the loop and continue execution 
until the last discontinuity is removed. 

Here we chose a slightly different approach, which can 
be implemented by a very minor change in the origi- 
nal formulation. Instead of picking an initial entrance 
spin and changing its state we leave its state unchanged. 
Without a magnetic field the loop would close immedi- 
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FIG. 4: A loop update in the presence of a transverse field. 
Note that there is no initial link-discontinuity, since the state 
of the initial spin in (a) is left unchanged. With a transverse 
field the first and last spins in the loop do not necessarily 
coincide. 



ately. With a magnetic field this possibility still exists, 
but the loop may also start with a flip vertex and there- 
after proceed in a "normal" fashion until another flip 
vertex is inserted, or removed, when the loop is closed 
according to the usual criterion (the exit spin and the 
following entrance spin are in the same state). We have 
compared high-precision QMC data for small systems 
with exact diagonalization, to ensure that detailed bal- 
ance is indeed satisfied, but we do not provide a strict 
proof here. 

Fig. ^ depicts an example of a loop move for the case 
of a non-zero transverse field. In Fig. ^(a) an initial en- 
trance spin is selected and left in its original state. An 
exit spin is selected and flipped (Fig. §(b)), resulting in 
a flip vertex. The next entrance spin can be seen in 
Fig. ^(c), with a corresponding exit spin in Fig. ^(d). 
Note that the exit state in Fig. ^(d) is unchanged, result- 
ing in a second flip vertex, leading to termination of the 
loop. The loop does not close on itself in this case, but 
starts and ends at two different spins. 

The transverse magnetization is particularly easy to 
calculate within the SSE formulation, since it simply is 
equal to the averase, number of flip vertices, iVflip, in the 
operator sequence,ll3 



(12) 



To conclude this section we note that only two changes 
need to be made to add a transverse field. First of all 
the initial entrance spin state should be left unchanged. 
Secondly, the exit probability is now also a function of 
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exit state, since there are, in general, several possible exit 
states for a given exit spin. Next we apply this algorithm 
to a two-dimensional ferromagnetic system. 



III. APPLICATION TO THIN MAGNETIC 
FILMS 

The simplest effective model for a thin magnetic film 
is a ferromagnetic Heiseriherg monolayer. However, the 
Mermin- Wagner theoremE£l tells us that this model can- 
not have a finite critical temperature in two dimensions. 
The continuous symmetry can be explicitly broken by 
an anisotropy in order to induce a finite critical tem- 
perature. Typical experimentally observed values of the 
local anisotropy in 3d transition metal films are about 
two orders of magnitude smaller than the exchange cou- 
pling, and one might therefore expect a very small criti- 
cal temperature. However, it turns out that the critical 
temperatuife,-e©Htains a logarithmic dependence on the 
anisotropy,Elt3EZI so that very small anisotropics induce 
a critical temperature of the order of the coupling con- 
stant. In this section we concentrate on results for the 
magnetization as a function of temperature. The model 
is given by Eq. where the double sum is over all 

nearest neighbours on a two-dimensional square lattice. 
Since the main purpose of the present paper is to show 
the feasibility of the quantum Monte Carlo method, and 
not to explore the whole parameter range of the model as 
defined in Eq. (pT|), we limit ourselves to some illustra- 
tive examples. From now on we will assume an isotropic 
exchange interaction (J*j = zero vertical mag- 

netic field {B^ — 0) and only second-order anisotropy 
[Ka^ — 0). First we discuss how the finite-size effects 
were treated. Thereafter we compare the QMC results 
with approximate theoretical approaches. Finally, we 
also show a case where a transverse field drives the mag- 
netization in the plane. 



A. Finite-size effects 

The first problem one encounters in a finite-size sys- 
tem with spin-inversion symmetry is that the magneti- 
zation should vanish because opposite spin orientations 
occur equally likely. One can circumvent this problem 
by calculating either the absolute value of the magneti- 
zation, or the magnetization squared (and take the root 
afterwards). Both approaches are equivalent in the ther- 
modynamic limit, and, since the absolute value of the 
magnetization is found to converge faster, we show only 
the former. The magnetization curves for a particular 
value of the anisotropy, K^j J = 0.01, and different sys- 
tem sizes are shown in Fig. ||. The finite-size effects in- 
crease closer to the Curie temperature Tc, as expected 
for a second-order phase transition. For low tempera- 
tures, however, one can see that it is possible to extract 
data which have converged (within statistical error) with 
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FIG. 5: Magnetization as a function of linear system size A'^ 
and temperature T/J for S = 1 and K2/J ~ 0.01. 
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FIG. 6: The Binder ratio Q as a function of linear system 
size A*' and temperature T/J for 5 = 1, and K2/ J = 1.5. 



respect to system size. In the next section we restrict 
ourselves to results for the magnetization that show no 
discernible finite-size effects. 

The critical temperature can-be determined directly. 
We have used the Binder ratioc3 to extract the value of 
Tc- By plotting the ratio of two moments of the magne- 
tization, 



Q 



(13) 



the finite-size effects around Tc should largely cancel, and 
curves for different system sizes are expected to intersect 
at Tc- In Fig. || we show an example for K2/ J = 1.5, 
where one can clearly see how data for different system 
sizes intersect at one point (Tc/J). 

In transition-metal thin films, the second-order 
anisotropy is beheved to be of the order K2/J — 0.01. 
The relatively small energy scale of the anisotropy in- 
duces large finite-size effects. In order to see the effects 



6 



O 





0^5 T^' 



-•QMC \ 

♦ QMCT^\ 
0(N) \ 

- SU(N) 

- RPA 
-- MF 



,/J 1.5 2 2.5 3 
T/J 



FIG. 7: The Binder ratio Q as a function of linear system 
size A'' temperature T/J for S = 1 and K2/J = 0.01. 



FIG. 8: Magnetization as a function of temperature T/J for 
S = 1 and K2/J = 0.01. 



of a small energy scale one needs to examine large sys- 
tems at low temperatures. In Fig. |^ we show the Binder 
ratio as a function of temperature and system size for 
K2/ J — 0.01. It seems that results for linear system sizes 
= 4, 8 and 16 have converged around Tc/ J = 1.05, but 
as the system size is increased strong corrections appear 
and push Tc/ J down to around T^/J = 0.9. This is a 
clear case where it would be dangerous to draw conclu- 
sions from a study of small system sizes. The results for 
system size = 256 have statistical errors which are 
slightly larger than the symbol size, while all other sta- 
tistical errors are much smaller. To determine Tc very 
accurately for such a small anisotropy one would need 
accurate data for even larger system sizes. However, the 
current precision is enough for the comparisons with re- 
sults of approximate theoretical methods that we present 
in the next section. 
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FIG. 9: Magnetization as a function of temperature T/J for 
5 = 1 and K2/J = 0.20. 



B. Magnetization 



We 



he QMC data with Schwinger Bo- 
and many-body Green's functions 
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calculations .ESEiJ Both methods include spin wave exci- 
tations approximately and represent significant improve- 
ments over simple mean field theories in which magnon 
excitations are neglected completely. Of these two meth- 
ods the Schwinger boson theories are numerically less de- 
manding and much easier to extend to arbitrary spin than 
the RPA approach .Ej 

In the Schwinger boson theory the Heisenberg model is 
mapped onto an equivalent bosonic system. This can be 
done by using the SU(2) symmetry in spin space of the 
Heisenberg model. The SU(2) model can then be gener- 
alized to an SU(7V) model, containing N bosons per site. 
In the limit TV — s- cx) mean-field theory becomes exact, 
and in this section we label mean- field SU(iV 00) re- 
sults by "SU(7V)" . Using the local equivalence between 



the SU(2) and 0(3) groups, this can be repeated for an 
0(3) model, and we label mean-field 0{N — > 00) results 
by "0(A^)". 

The many-body Green's function calculations for the 
magnetization are done by a procedure where one works 
at the second level of the hieraphy of equations of mo- 
tion for the Green's functions.tll This allows an exact 
treatment of the terms stemming from the single-ion 
anisotropy, whereas thCpexchange interaction terms are 
treated by the TyablicoMlj, or random-phase approxima- 
tion (RPA), and from now on Green's function results 
will be labeled by RPA. This procedure is ran improve- 
ment over the Anderson-Callen decouplingll3, in which 
the single-ion anisotropy terms are decoupled at the level 
of the lowest-order equation of motion, which is a good 
approximation only for small anisotropics. til The results 
in this subsection are calculated using the exact treat- 
ment of the single-ion anisotropy. Owing to problems 
with numerical stability, the calculations for the reorien- 
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FIG. 10: Magnetization as a function of temperature T/J 
for S' = 1 and K2/J = 1.50. 



tation of the magnetization discussed in the next subsec- 
tion are done with the Anderson-Callen decouphng. 

We compare magnetization curves for three values 
of the anisotropy covering three orders of magnitude, 
^2 /■/ = 0. 01, 0.2 and 1.5. The smallest value is of great- 
est relevance for experiments, but it is of interest to see 
how well the analytic methods work outside this region 
as well. In Fig. | we see that, for K2/J = 0.01 the RPA 
calculations give a rather accurate magnetization curve, 
while the 0{N) and SU(iV) theories give a low and high 
estimate, respectively. In the low-temperature limit both 
the SU(iV) theory and the RPA calculation recover the 
correct spin- wave result. Note also by how much a sim- 
ple mean-field (MF) theory overestimates the magnetiza- 
tion {T^^ ~ 2.7 TQ^c !). As is well known, mean-field 
theory also totally fails at low temperatures (exponen- 
tial instead of power law behaviour), due to the neglect 
of spin waves. In Fig. ^ we see that for K2/J = 0.2 
the RPA still gives very accurate values at low tempera- 
tures, while overestimating the magnetization at higher 
temperatures. The 0(7V) theory happens to give a good 
estimate of Tc, while the SU(7V) again overestimates the 
magnetization. In Fig. ^ the QMC points are connected 
by straight lines, except for the line between and the 
highest temperature below Tc- This line is a fit to the 
Ising-like critical behavior 



(14) 



The good fit indicates that the QMC data have come 
close enough to Tc for the critical behavior to set in. For 
the largest anisotropy, K2/J ~ 1.5 (Fig. |l^), the RPA 
calculation again gives the most accurate estimate, but 
all approximate curves result in too high a magnetization. 
For high temperatures the SU(A^) theory yields a lar ger 
magnetization than simple mean- field theory. In Fig. |lO| 
we have also included the quantum Monte Carlo data 
raised to the eighth power. The resulting straight line 
is further evidence of how close to the critical point the 
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FIG. 11: In-plane and out-of-plane magnetization as a func- 
tion of transverse field at temperature T/J = 1 for S = 2, 
K2/J = 0.01. 



simulation has come. 

The large deviations of the approximate theoretical 
treatments from QMC, in particular close to Tc and for 
large anisotropics, indicate that neither the Green's func- 
tion approach nor the Schwinger boson method treats 
spin wave interactions in a satisfactory way. It is interest- 
ing to note that that RPA and SU(iV) give a (mean- field) 
exponent of 1/2, whereas O(A^) gives 1/3. 

In this subsection we have compared QMC data with 
approximate theoretical results for the magnetization 
curve. Except for the value of Tc, the QMC data shown 
above have converged in system size within error bars 
that are not discernible in the figures. The largest sys- 
tem size used is 256 x 256 spins. Error bars are shown 
for the Binder estimate of Tc. As can be seen the er- 
ror increases with decreasing anisotropy, which is to be 
expected. 



C. Reorientation 

We consider two examples for the reorientation of the 
magnetization in a transverse field. For an anisotropic 
spin-2 model we have calculated the vertical and trans- 
verse components of the magnetization, and M^., as 
a function of a transverse magnetic field at a fixed tem- 
perature. 

The anisotropy favors an out-of-plane magnetization, 
while the transverse field wants the magnetization to be 
in the plane. This competition results in quite interest- 
ing phase diagrams, where the order of the transition in 
general depends on the order of the anisotropyJl3 Here 
we concentrate on two examples, to demonstrate the ap- 
plicability of the QMC approach. 

In Fig. O we show the magnetization curves for a small 
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FIG. 12: In-plane and out-of-plane magnetization as a lunc- 
tion of transverse field at temperature T/J = 1 for S = 2, 
A'a/J = 0.2. 

anisotropy, Kij J = 0.01, at a temperature T/J = 1. Pri- 
marily we show QMC results for system size 128 x 128, 
but in addition we show results for 64 x 64 in the region 
where they differ from each other. For small and inter- 
mediate fields the results have converged, while we again 
see increasing finite-size effects closer to the critical field. 
The trend is, however, clear. The transverse magnetiza- 
tion does appear to increase linearly with the transverse 
field, up to the highest fields for which the results have 
converged. 

The SU(A^) and O(A^) calculations both show a linear 
dependence on the transverse magnetization, M^^ similar 
to QMC, whereas the vertical component is over- or un- 
derestimated, respectively. The SU(A^) results for are 
good at weak field, though. The mean-field results for 
are also linear but mean-field theory consistently overes- 
timates Mz ■ The RPA calculation with the approximate 
decoupling of the anisotropy termsEJ follows the QMC 
curves at small values of the transverse field, but the re- 
orientation occurs at a considerably smaller critical field 
than in QMC and in the other approximations. Also, in 
the RPA the magnetization is not a unique function of 
field close to the transition. The same system has been 
studied for a larger anisotropy, K2/J = 0.20, in Fig. [l^ . 
Here the SU(A^) and 0(iV) solutions are mean- field- like 
and agree very well with the QMC solution, except close 
to the transition, whereas the RPA behaves as in Fig. |l^. 
Close to the transition the SU(A^), O(A^), and mean- field 
solutions overestimate the vertical magnetization. The 



QMC results display a small deviation from the linear 
increase in the in-plane magnetization, which is not re- 
produced by the Schwinger boson or mean-field methods. 
IV. CONCLUSION 

This work shows the feasibility of using large-scale 
QMC calculations to examine microscopic thin film mod- 
els. The QMC approach can be used both as a san- 
ity check on approximate theoretical treatments and as 
a method in its own right. Our results indicate that 
in the absence of a transverse field- RPA with an exact 
treatment of the anisotropy termstil appears to be more 
accurate than the Schwinger boson calculation. Both 
methods, however, represent an improvement over sim- 
ple mean-field theory for small anisotropics but QMC 
reveals weaknesses in the approximate theories at large 
anisotropics and close to the Curie temperature. In 
the presence of a transverse field the results obtained 
from Schwinger boson methods are quite mean-field-like, 
which turns out to be appropriate for most of the field 
range. The reason that the RPA is worse in this case 
is probably due to the Andepaon-Callen decoupling of 
the anisotropy terms in RPA .113 This_was necessary be- 
cause the more accurate treatmentliil led to numerical 
difhculties when applied to the reorientation problem. 
To make more definite conclusions about the merits of 
different approximate methods, we would have to inves- 
tigate a much larger parameter space, including higher 
anisotropics, general spin, temperature, and transverse 
field and the extension to several layers. We leave such 
an investigation as work for the future. 

It is possible to extend the QMC calculations to include 
several layers with arbitrary inter-layer coupling, as well 
as, for example, anisotropics in the exchange coupling. 
Unfortunately, the dipole interaction introduces frustra- 
tion and thereby the sign problem, and is therefore cur- 
rently out of reach for QMC studies. On the other hand, 
long-range ferromagnetic couplings are not a problem. 
The implementation of the dipole coupling also leads to 
problems in the Schwinger bosMi-t|heory but is possible in 
a Green's function description.Qty A very exciting recent 
develoBHient in QMC is the introduction of directed loop 
movcsJHJ which, according to our initial calculations, can 
reduce the autocorrelation time by one order of magni- 
tude in spin-1 systems, and therefore make it possible to 
reach substantially larger system sizes. However, it is not 
clear that the directed loops are as easy to implement for 
a general model as the method used in this work. 

We are grateful to A. Sandvik for fruitful discus- 
sions. P.H. acknowledges support by the Swedish Re- 
search Council and Ella och Georg Ehrnrooths stiftelse. 
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